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(54) Method and apparatus for eliminating background interference signals from multichannel 
detector arrays 



(57) A method and apparatus is described for mak- 
ing multichannel measurements of weak signals in noisy 
environments. The method is based on an adaptive 
compensation technique in which the large interfering 
background signal is first recorded by using the appa- 
ratus itself. By a statistical analysis of this multichannel 
measurement the independent components of the inter- 
ference are determined. The apparatus is provided with 



compensator elements which are properly coupled to in- 
dividual sensors for both to collect and distribute infor- 
mation on the interference during the measurement. 
The proper couplings are derived from the interference 
measurement mentioned above. In this way the output 
of the device can be made immune to the large ampli- 
tude interference components present in its environ- 
ment. This reduces the dynamic range requirement for 
the data transport and storage systems. 
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Description 

FIELD OF THE INVENTION 

5 [0001] The present invention relates to multichannel measuring. More particularly, the present invention relates to a 
novel and improved method for collecting multichannel signals comprising of the signal of interest and of superposed 
background interference contributions which may be much larger than the signal of interest. 

DESCRIPTION OF THE RELATED ART 

10 

[0002] Performing many simultaneous measurements on a subject, i.e. multichannel detection, is sometimes essen- 
tial in order to obtain sufficient information on the issue under examination. We consider, in particular, the detection of 
biomagnetic fields associated with the function of human brain or heart. Modern magnetometers for this purpose com- 
prise about 100 channels to enable accurate localization of neuro- or cardiography sources. Biomagnetic fields are 
15 very weak in comparison to the background magnetic fields in the surroundings, so that the problem of resolving the 
real signal from environmental interference is technically very challenging (M. Hamalainen et. al., "Magnetoencepha- 
lography - theory, instrumentation, and applications to noninvasive studies of the working human brain", Rev. Mod. 
Phys. vol. 65, no 2 April 1993.). 

[0003] Prior art of protecting very sensitive instruments against external interference includes basically five methods: 
20 1) use of passive shielding elements surrounding the instrument (magnetically shielded room in the biomagnetic ap- 
plication), 2) use of active elements canceling the interfering environmental signal (large scale compensation coils in 
magnetic measurements), 3) reducing the relative sensitivity of the sensors to typical background signals (use of gradi- 
ometers instead of magnetometers), 4) use of additional sensors to estimate the background interference in order to 
separate it out from the signals, and 5) numerical processing of the multichannel data to separate true signal from 
25 external interference. 

[0004] In method 1), when applied to biomagnetic measurements, the instrument is placed inside a shielding room 
having walls made of high permeability metal alloy (mu-metal). In the low frequency range, relevant to biomagnetic 
signals, the shielding factor of such a room is limited to about 100 - 1000 by reasonable amount and finite permeability 
of mu-metal. At high frequencies the shielding may be improved by adding layers of highly cbnducting material, such 
30 as aluminium (V. O. Kelha et al., "Design, Construction, and Performance of a Large-Volume Magnetic Shield", IEEE 
Trans, on Magnetics, vol. MAG-18, no 1, Jan. 1982.). 

[0005] In studies of human subjects, possibly patients in a hospital, the magnetically shielded room has to be relatively 
large, leading to a heavy and expensive construction. Sufficient shielding requires multilayer structure with total wall 
thickness of about 0.6 m. Thus, the outer dimensions of the room must be on the order of 4 m x 5 m x 3.5 m to provide 
35 enough space for the instrument and comfortable conditions for the patient on a bed, and possibly for medical personnel 
taking care of the patient. Especially, the need of 3.5 m in height (two floors) is inconvenient in a typical hospital 
environment. 

[0006] Method 2), when large compensation coils are used (EP 0 514 027, M.Kazutake et al. "Magnetic noise re- 
ducing device for a squid magnetometer") resembles the passive shielding with high permeability material. The shield- 

40 ing current, which in mu-metal is generated as a natural response to an exposure to magnetic field, is now generated 
artificially in a control system and driven into coils with dimensions comparable to those of a typical shielded room. As 
a realization of such a system, three orthogonal Helmholtz pairs may be used. The external field to be eliminated is 
measured outside the coil system by field sensing elements, such as fluxgates, whose output is converted by a proper 
control system into electrical currents fed into the compensation coils. This kind of active shielding is far lighter and 

45 less expensive than a typical passive shield. It also performs best at low frequencies, where passive shielding of 
magnetic fields is most difficult. 

[0007] The major disadvantage of method 2) is the very restricted geometry of the shielding currents. In practice, a 
compensation coil system can reject the field of distant sources only, which produce nearly uniform field at the site of 
the instrument. It may also be difficult to find the optimal positions for the field sensing elements, and if the environmental 

50 conditions change, the system may have to be readjusted. 

[0008] Method 3), regarding the biomagnetic application, is based on the fact that the gradients of a magnetic field 
decrease more rapidly as a function of the distance from the source than the field itself. Therefore, the signal to back- 
ground ratio is increased by measuring the difference of magnetic flux between two adjacent locations instead of the 
flux itself: the signal arising from the nearby object of study (e.g. a brain) is enhanced in comparison to the disturbance 

55 signal from an interfering source further away. 

[0009] In principle, method 3) provides total immunity against uniform interference fields. In practice, however, the 
balance of best gradiometers is limited to at best 1/1000 because of technical difficulties in controlling the geometry 
of the sensors. In addition, the interfering fields are never strictly uniform. If the disturbing source is located a distance 
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1 away (typically 1 - 10 m) and the baselength of the gradiometer is h (typically 0.01 - 0.1 m), the background signal 
of the sensor is damped roughly by a factor of h/l compared to a magnetometer with the same loop size. 
[0010] The most severe drawback of method 3) is that it rejects part of the signal arising from the object of study as 
well. This is especially unfavorable when the biomagnetic field is nearly uniform on the length scale of the sensor. This 
5 is to some extent the case in cardiac studies and when a neuromagnetic source is located deep below the scull. For 
this reason, magnetometers would be preferred instead of gradiometers in many biomagnetic measurements (M. 
Hamalainen et. al., "Magnetoencephalography - theory, instrumentation, and applications to noninvasive studies of the 
working human brain", Rev. Mod. Phys. vol. 65, no 2 April 1993). 

[0011] In method 4) ( US 5,187,436 A, J. A. Mallick "Noise cancellation method in a biomagnetic measurement 
10 system using an extrapolated reference measurement", and US 5,020,538, N. H. Morgan et al., "Low Noise Magnet- 
oecephalogram system and method", and DE 41 31 947, G. M. Daalmans, "Mehrkanalige SQUID- Detektionseinrichtung 
mit Storfeldunterdrockung", and DE 430451 6, K. Abraham-Fuchs, "Verfahren zum Bestimmen einer Characteristischen 
Feldverteilung einer ortsfesten Storquelle", and WO 93/17616, K. Abraham-Fuchs, "Disturbances suppression process 
during position and/or direction finding of an electrophysiological activity", and EP 0481 211 , R. H. Koch, "Gradiometer 
15 having a magnetometer which cancels background magnetic field from other magnetometers", and US 5,657,756, J. 
Vrba et al., "Method and systems for obtaining higher order gradiometer measurements with lower order gradiometers") 
the apparatus is equipped with additional background sensors, which are so arranged that they do not receive any 
substantial input from the object of study. They are usually placed further away from the actual sensor array. From the 
signals of these sensors an estimate of the interfering background field is calculated - for example up to the desired 
20 order in the Taylor expansion of the field - and then properly extrapolated and subtracted from the signals of the actual 
measuring channels. 

[0012] The relatively large distance between the background sensors and the actual sensors and the inaccuracy in 
the calibration and relative location and orientation of the sensors are the main drawback of this method, because 
these factors limit the degree of achievable compensation. Especially, correct interpretation and use of the background 
25 sensor outputs is practically impossible, for example, if the background signal arises from an unknown vibration mode 
of the instrument in an unknown remanence field distribution. 

[001 3] In method 5) the signals collected by a multichannel device during a measurement are first stored on a memory 
device. After the measurement, the data are processed with a numerical template or projection method to separate 
out the contributions of the interesting biomagnetic sources from the disturbing interference fields (WO 94/12100, R. 
30 llmoniemi, "Method and apparatus for separating the different components of evoked response and spontaneous ac- 
tivity brain signals as well as of signals measured from the heart", and WO 93/1 761 6, K. Abraham-Fuchs, "Disturbances 
suppression process during position and/or direction finding of an electrophysiological activity", and US 4,977,896, S. 
E. Robinson et al., "Analysis of biological signals using data from arrays of sensors"). 

[0014] Method 5) relies on the multichannel aspect of the measuring device: only by collecting data from many 
35 sensors simultaneously can the background interference be separated from the true signal due to their characteristically 
different distribution over the entirety of channels. In a neuromagnetic measurement, for example, the sensors should 
cover the whole head. 

[0015] When used with magnetometers method 5) requires very large dynamic range for every channel of the data 
collection system, since the actual biomagnetic signal can be contaminated by a background signal several orders of 
40 magnitude larger. One would effectively have to subtract large but nearly equal numbers from each other to reveal the 
differences representing the actual biomagnetic activity. 

[0016] In practice, the strength of biomagnetic fields is 6 - 8 orders of magnitude weaker than the unshielded back- 
ground fields in a typical environment (M. Hamalainen et; al., "Magnetoencephalography - theory, instrumentation, and 
applications to noninvasive studies of the working human brain", Rev. Mod. Phys. vol. 65, no 2 April 1993). Therefore, 

45 at least two of the above methods have to be combined to achieve a tolerable signal to background ratio. 

[0017] In the present invention a multichannel sensor device is made immune to environmental interference by cross 
coupling the channels in such a way that there is no output in response to the interference. No extra compensation or 
reference channels are necessary. In practice it turns out that the shielding efficiency of the present method is propor- 
tional to the number of cross coupled channels and therefore the present method where all or most of the signal 

50 channels participate the compensation is supenor to prior art methods which utilize a smaller number of separate 
compensation channels (e.g. US 5,657,756, J. Vrba et al., "Method and systems for obtaining higher order gradiometer 
measurements with lower order gradiometers"). 

[0018] Ordinarily, when N channels are operated in parallel, the output of each channel depends on the input of its 
own sensor only. This can be described by a diagonal NxN matrix C, by which the N-dimensional output vector U is 
55 obtained as 



U=Cu 
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for a given N-dimensional input vector u. The element 0 {] of the coupling matrix is the gain, or the calibration constant, 
of the respective channel /. 

[0019] The present compensation method is described by a non-diagonal matrix C, whose off-diagonal components 
5 represent the cross couplings between the channels. This matrix is constructed so that it maps to null vector all the 
input vectors interpreted as interference; the required linear mapping C in N-dimensional signal space has n-dimen- 
sional null-space, where n is the number of independent interference vectors, or field distributions, spanning the sub- 
space called interference signal space. In practical applications n«N. 

[0020] The prior art invention WO 94/12100 (R. Ilmoniemi, "Method and apparatus for separating the different com- 
10 ponents of evoked response and spontaneous activity brain signals as well as of signals measured from the heart") 
describes a signal space method where different biomagnetic responses are separated from each other and from 
interference originating from uninteresting sources by applying signal space projection methods to collected data. 
Typical interference signals - especially in the case of magnetometers - may be by factor 10000 or 1000000 larger 
than the signals of interest. Therefore the aforementioned prior art method would require data collection and storage 
15 with too much extra dynamic range to be practical. 

[0021] The required cross-coupling strengths for the sensor network in the present method are determined from a 
measurement of the interference seen in the absence of the cross couplings. For the determination of the cross coupling 
strengths no detailed information on the location, orientation, or calibration of individual channels or their relative lo- 
cations and orientations is required. Full compensation of an A/-channel system can be accomplished by 2Nn cross 
20 couplings, when the couplings are realized by negative feedback. In practice even a lower number of couplings may 
be sufficient. 

[0022] When performing a measurement with the compensated system the magnitudes of the background interfer- 
ence signal components in the interference signal space are recorded together with the compensated signals and, if 
required, the uncompensated signals can be recovered from this information with a linear transformation. 

25 [0023] The use of the present method in the biomagnetic application is in a way analogous to using a shielded room: 
The sensors of the multichannel magnetometer detect the interference due to the field left in the magnetically shielded 
room and n "shielding currents" are constructed from this information and then delivered, properly weighted, to the 
individual channels in form of negative feedback. This negative feedback is superimposed on the ordinary negative 
feedback used to drive the magnetometer channels in the flux locked loop. The same feedback coils used for the flux 

30 locking negative feedback can be used to feed in the shielding currents as well. 

[0024] As to prior art methods 1) and 3) the invention effectively improves magnetic shielding and thus enables to 
use sensitive magnetometer sensors instead of gradiometers in a standard shielded room. The present method for 
eliminating the interference is adaptive to the conditions present at a particular site, since the cross couplings are 
chosen to cancel the interference measured by the very sensor array itself. The effective shielding factor so achieved 

35 js comparable to that of the best balanced gradiometers of the date. 

[0025] As to prior art method 2) the present invention offers more flexible adaptive shielding. The N negative feedback 
coils of the individual channels replace the small number of large, fixed geometry external compensation coils of method 
2). In our method the compensation currents have no interaction with the walls of the shielding room. Also, any vibration 
of the magnetometer array in the remanence field is impossible to handle with a set of external compensation coils but 

40 in the present method it is simply an extra dimension in the interference signal space. 

[0026] As to prior art method 4) the present invention offers several advantages: No extra compensation or reference 
sensors are necessary. For successful compensation there is no need to accurately calibrate or balance sensors or to 
make any sensors parallel or orthogonal to each other. Neither is the compensation limited to any order in the Taylor 
expansion of the interfering magnetic field. In the present method the interference is simply compensated up to any 

45 order necessary; the degree of compensation achieved depends only on how accurately the devices used for setting 
the cross coupling strengths can be set. The degree of achievable compensation is also proportional to the number of 
channels participating in the cross coupled network. In the present method this number can be freely chosen and can 
be increased up to the total number of channels in the device (N) instead of the relatively limited number of separate 
compensation channels used in the prior art compensation schemes. A standard compensation obtained by adding or 

50 subtracting a reference signal increases the noise in the signals. Such increase of noise is absent from the present 
method because the linear mapping C is a projection. 

[0027] A reduction in the number of actual measurement channels takes place in method 4) when out of the N sensors 
of the system n are chosen permanently to be compensation sensors and moved further away from the source of the 
actual signals. The present method is more flexible at this point because the number of compensated interference 
55 modes n can be chosen according to the needs dictated by the environment and the quality of the shielded room. 
[0028] Also, in method 4) the quality of compensation may suffer if any one of the n compensators does not work 
properly. In the present method this problem is absent because the individual channels are equal and a malfunctioning 
channel can simply be disconnected from the feedback network (and discarded from the data). 
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[0029] As to prior art method 5) the advantage is that no sensors or data collection devices with excess dynamic 
range are needed. The "shielding currents" created in the cross coupled sensor network are distributed to balance all 
the N sensors against the interference. 

SUMMARY OF THE INVENTION 

[0030] The object of the present invention is to eliminate the problems and the disadvantages described above. 
[0031] A specific object of the present invention is to disclose a completely new type of method and device for elim- 
inating background interference from multichannel signal measurement. 

[0032] The approach of the present invention can be used in conjunction with measurements by any multichannel 
device, susceptible to environmental interference. The preferred embodiment is directed toward an application for 
biomagnetic measurements: magnetoencephalography (MEG) and magnetocardiography (MCG). The method, how- 
ever, is more generally applicable, as long as the environmental background signal has sufficiently different charac- 
teristics from the signal of interest. This condition can usually be fulfilled by a properly arranged sensor array with a 
sufficiently large number of channels. For example, in neuromagnetic measurements this is accomplished by a whole 
head coverage with about 100 sensors. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0033] In the following, the invention and its advantages are explained in more detail with reference to the attached 
drawings, in which 

Fig. 1 presents a schematic illustration of a simplified sensor arrangement in a multichannel measurement, 

Fig. 2 presents a schematic illustration of a signal processing in a simple cross coupled multichannel device, and 

Fig. 3 presents a schematic illustration of a two measuring channels and a compensating channel with a negative 
feedback network realized by means of an operational amplifier 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

[0034] The present invention is most easily described by means of a highly oversimplified example, which, however, 
is straightforwardly generalized to apply in actual practical implementations. Therefore, we consider a device with only 
two sensors; let them be magnetometer loops for clear visual interpretation, see Fig. 1 . Let these loops, (1 and 2), be 
in the same plane, let the interference be a uniform field in arbitrary direction, and let the signal of interest arise from 
a current-dipole like source (3) causing the magnetic fluxes <)) 1 and <j> 2 through the two loops. 
[0035] The input vector u is the sum of the uniform background flux % penetrating both loops plus the contributions 
of the source of interest, <|> 1 and (<t> 2 . When the two channels are read separately, i.e. the system is diagonal, the matrix 
C is a 2 x 2 unit matrix multiplied by an overall calibration constant C 0 . The input vector u =(<|>o + <|>i, (^o" 1 "^) results in 
the output vector U= Cu = C^o*^, <t>o + <l>2)* ^ the background is very large compared to the actual signal of interest, 
the output is dominated by the interference and the actual signal from the dipole can not be resolved unless the sensors 
have a very large dynamic range. 
[0036] However, if the coupling matrix is modified to 



we obtain U = (C^2)(^-^ 2 * ^i +< t>2) f° r tfie same input vector as above. Now the device is made completely insensitive 
to the large uniform interference field and the response is only due to the dipole source. The price paid for this con- 
venience is that the two signals have become linearly dependent (equal but of opposite sign). For large sensor arrays 
this is not a problem, since compensation of n independent interference components from the signals of N channels 
reduces the number of independent signals to N- n. When N » n, there is no essential loss of signal dimensionality. 
[0037] The non-diagonal couplings mean that the input of each channel contributes to not only the output of the 
channel itself but also to that of the other channels. The parallel, independent sensors of a conventional device are 
thus transformed into an integrated sensor network. This is illustrated in Fig. 2. 




i -fi 
-i ij ' 
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[0038] The example above can be generalized to a device of N magnetometer loops, all laying still in the same plane 
for simplicity. A signal in response to a uniform background field is expelled by using an N x N coupling matrix 



10 



15 



20 



(N-l -1 
-1 N-l 

I -1 -I 



-1> 
-1 

N-l) 



■ (1) 



Here the output signal of each channel consists of the diagonal input, weighted by (N - 1)/A/, and of the (N- 1) contri- 
butions with weights -1//V from all the other channels. 

[0039] In this scheme a multichannel device with N sensors is compensated against an external disturbance by 
making N 2 cross couplings. It would become exceedingly complex and expensive to realize such an amount of cross 
couplings for systems consisting of about 100 channels. This complication is avoided when the cross couplings are 
made as described in the present invention. 

[0040] When the two channels in the example above are operated in a negative feedback loop, they may be cross 
coupled to receive feedback also from each other 



25 



[0041] As in ordinary negative feedback the forward gain is assumed to be virtually infinite i.e., G 1t G 2 » Kyi /, j = 
{1 , 2}, and we end up with a matrix relation u= -kU, where the matrix k has the elements 1/Ky. The output vector is U- 
30 -k- 1 w , i.e. the linear mapping C above is defined by the matrix 4r 1 . 

[0042] The uniform background interference is canceled from the two channels when the non-diagonal couplings 
are chosen to give 



35 i ( 1 -ft 

k -vU J ■ 

[0043] This matrix, however, is singular. Consequently, no realizable negative feedback couplings exist, because 

40 they are determined by the inverse of the coupling matrix. 

[0044] This problem can be overcome by the present method. A "virtual channel" or element is added to the system. 
This channel has no actual sensing element; it simply acts as a stage collecting and distributing the feedback between 
the real channels and appears as an additional component in the signal vector. In the following, this channel will be 
called the compensator channel. Its output (U c in Fig. 3) is actually proportional to the projection of the input vector on 

45 the interference signal space and can be recorded and stored along with the compensated signals of the proper chan- 
nels. Within this scheme, the proper feedback matrix k in the above case of two magnetometers is 



50 




(0 1 






1 1 


0 






a o 


i, 



55 [0045] The first component of the signal vector corresponds to the compensator channel. With the corresponding 
coupling matrix 
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one can verily the desired result. The background input vector is now u = (0,1 .1 ) <t> 0 ( the compensator channel gets no 
input except from the other channels). The interference is not seen in the output of the actual signal channels; it is 
displayed in the compensator channel only: U = -k 1 (0,1,1)<t> 0 = -/^(I.O.O)^. 

[0046] The advantage of the cross-coupling scheme realized by using the compensator channel (or element) be- 
comes apparent when the number of channels is increased. In the case of N magnetometer loops in a plane, a uniform 
background is expelled from the measuring channels with the scarce ) x (A/+1) feedback matrix 



k = l/K n 
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(2) 



Couplings are needed only to and from the compensator channel. Instead of the N(N- 1) = 9900 couplings between 
100 channels it suffices now to make 2N = 200 cross couplings. 

[0047] The compensator channel can be realized, for example, by a simple adding operation amplifier, as illustrated 
in Fig. 3. The output of each real channel is connected to the input of the said amplifier (4), the output of which is 
connected again to the feedback line of each channel. Each channel's own negative feedback and the cross coupling 
feedback from the compensator channel are added up, for example by another operation amplifier (5 and 6), and 
coupled to the sensor elements (1 and 2) via the respective feedback coils. In this realization the elements of the 
feedback matrix, i.e. the cross-coupling constants, are determined by the input resistors of the operation amplifiers. 
The "forward couplings", i.e. the off-diagonal elements on the first row of the feedback matrix, are scaled by the inverse 
of the resistors R 12 and R 13 ; the "reverse couplings", i.e. the off-diagonal elements on the first column of the feedback 
matrix, are scaled by the inverse of the resistors R 2 1 aR d ^31- The overall scaling of the matrix elements is determined 
by the feedback resistors R f . They are preferably identical for each channel, but even if they are not, the adaptive 
nature of the present background elimination scheme will take care of that also. 

[0048] In a system with digital feedback loop the cross couplings can be realized numerically by a signal processor. 
The analog output of a magnetometer channel is converted into digital form and then combined, properly weighted in 
the processor, with the outputs of the other channels. 

[0049] The possibility to use information gathered with a large number of magnetometers (A/) — instead of the infor- 
mation collected with just a small number of reference channels — is advantageous for the present interference com- 
pensation scheme. This becomes obvious when looking at the frequency dependent responses. Assuming that the 
bandwidths of the cross-coupling operational amplifiers ((1 ) in Fig. 3) are given by G c = 1/io>x cl and the transfer function 
of the diagonal feedback loop is given by G/K = -l/ion, the output of the actual channels in response to an input signal 



50 



(7j = ±K/(1+io)T)y s , 



55 



and the output in response to the interference 

U. = -K/(1-icot+A//(icox c ))u 0 . (3). 



[0050] The response of actual channel to signal is unaffected by the compensator whereas the response to low 
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frequency interference (co « r 1 , x c " 1 ) is reduced by an effective shielding factor S = u^U^ = NI(Km c ) which is proportional 
to N. A large number of channels involved in the compensated sensor network implies high rate of compensation. 
[0051] According to Equation (3) the compensation system is least effective at co = (A//(xx c )) 1/2 where the shielding 
factor achieves a minimum value S min = (1 + 4/Vx/x c ) 1/2 . The minimum is pushed beyond the measurement band (w~1/x) 
5 by a factor (A/x/i c ) 1/2 , and even at the minimum the shielding factor is of the order (4/Vr/x c ) 1/2 . Thus, from this point of 
view a large number of channels in the compensation network is essential. 

[0052] The conclusions above on the time dependent response of the compensated system remain the same even 
if the channels are unidentical - only the analytic formulae become more complex. 

[0053] The generalization of the above compensation scheme to an arbitrary sensor arrangement and to any spatial 

10 form of interfering field is straightforward. Only the couplings, of strength unity above, are changed to values between 
-1 and +1 which — in the case of magnetometers — describe the projections of the magnetometer loop area along 
the direction of the interference field. Also, in the general case the number of compensated interfer-ence modes, i.e. 
the dimensionality of the interference signal space, must be increased beyond one used in the examples above. When 
the sensors of the array have arbitrary orientation — instead of being coplanar — even the compensation of an arbitrary 

15 uniform interfering field requires compensation of three components of the magnetic field (n=3). 

[0054] In a practical application, that is in a real shielded room, it is, however, not necessary to compensate the 
multichannel system against any ideal field distribution - like the uniform field, or the five independent first derivatives 
of the field etc. (see prior art; US 5,657,756, J. Vrba et al., "Method and systems for obtaining higher order gradiometer 
measurements with lower order gradiometers"). The essential need is to compensate against those field distributions 

20 observed to arise in the shielded room in response to the typical magnetizing events in the neighbourhood of the very 
room in question. If the room is located reasonably far from the major interference sources (streets, railways, elevators) 
the response of the room consists of three independent interference field distributions corresponding to the three com- 
ponents of the approximately uniform magnetizing field outside the room. These field distributions are practically never 
uniform inside the shielded room. 

25 [0055] The actual number of compensator channels is chosen according to the number of interference modes that 
require compensation. If the external field magnetizing the shielding room is essentially in one direction even one 
compensator may be enough. On the other hand if the walls of the shielded room or the magnetometer are vibrating, 
or if there are interfering signals arising from sources inside the shielded room (cardiac signal of the subject, or signals 
from electrical devices) the compensation of these interference components is done by adding one compensator chan- 

30 nel per interference mode. 

[0056] The method for determination of the cross-coupling weights, in a specific environment is an important aspect 
of the present invention. This is done by the sensor array itself without the cross couplings. No preknowledge of the 
nature of the interference is needed. The compensation is tuned to work in the conditions of the particular installation 
of the multichannel device. 

35 [0057] The interference signal is recorded in an empty shielded room with the door closed for typically a couple of 
minutes and a statistical analysis such as principal component analysis (PCA, see for example S. Haykin, Neural 
Networks, Macmillan College Publishing Company, pp. 363-394, 1994) or independent component analysis (ICA, see 
e.g. A. Hyvarinen and E. Oja, "A fast fixed-point algorithm for independent component analysis", Neural Computation, 
vol. 9, pp. 1483-1492, 1997) is made on this data. The dominant principal components of this empty room recording 

*o give a set of orthogonal basis vectors of the interference signal space. Prior to setting the cross couplings a numerical 
signal space projection (SSP, WO 94/12100, R. Ilmoniemi) method can be applied to verify the number of interference 
components necessary to compensate, so that the interference contribution in the output signal is acceptably small. 
[0058] In the vicinity of the shielded room there usually exists a dominant source of interference causing a magnet- 
izing field essentially in one direction. Therefore, the dominant interference mode found in the empty room signal PCA 

45 may well be an order of magnitude stronger than the next, weaker modes still to be compensated. In this situation, it 
may be practical to compensate some linear combinations of these principal components in order to load the -com- 
pensator channels more evenly and to gain more dynamic range. Two modes, c 1 and c 2 , can be combined to (Cj + C2) 
Nl and (c 1 - C2W2, for example, which retains their orthonormality. 

[0059] In practice, the adjustment of the coupling constants in a large scale system with about 100 channels must 
50 be computer controlled. Without direct readout of the set cross-coupling values, there must be means for checking the 
status of these couplings. This is necessary both for the success of the interference rejection and for the correct inter- 
pretation and analysis of the signals collected with the compensated device. After the cross-couplings have been set 
the actual coupling matrix can be verified column by column by exciting each channel at a time (including compensator 
channels). This corresponds to an input vector with just one element different from zero. In an ordinary uncompensated 
55 magnetometer array the hardware necessary for exciting one channel at a time already exists for tuning purposes. A 
standard way is to use the feedback coil of the channel for this purpose. Due to the cross couplings in a compensated 
system, the excitation of one channel alone generates an output signal to the other channels as well. The output vector 
produced by exciting i:th channel is directly proportional to the i:th column of the cross-coupling matrix. 
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[0060] If the actual cross-coupling matrix K which has been determined with the above method is found to differ from 
the desired cross-coupling matrix k it is corrected by making the element by element corrections Ay = (fcjj-[K 1 ],)/ [K _1 ]y, 
where }-{n^ , n+A/} for /= {1 , n) and j = {1 , n} for /= {n+1 , n+N}\ [K" 1 ]y denotes the ijith element of the inverse matrix of K. 
[0061] This procedure can be repeated until the compensation result is satisfactory. As a rule of thumb, to compensate 
5 an interference signal by an effective shielding factor S, the cross-couplings between the channels have to be set with 
a relative precision of S* 1 . 

[0062] Once the cross-couplings have been set, there is also a quicker way to check them at any time against acci- 
dental changes. A practically sufficient brief check is to measure the n x n submatrix by exciting and reading the 
compensator channels only. This submatrix must be a diagonal matrix A/ -1 1, when the coupling vectors n 4 are scaled 
10 as above, and the compensated interference vectors are orthonormal linear combinations of the principal components 
of the empty room recording. An unintended change in any feedback connection results in appearance of a non-zero 
off-diagonal element in the n x n submatrix. 

[0063] To ensure unbiased analysis of data recorded with the multichannel device such a bnef or full check of the 
cross couplings can be routinely made in the beginning or end of the measurement. 
15 [0064] The cross coupling strengths used to reject an interference component in the present method can be chosen 
in many different ways: It is not necessary to include all the channels in the compensation network. One can leave out 
the "forward" couplings (resistors R 12 , ft 13 in Fig. 3) of any subset of channels. This allows one to prevent occasional 
noisy channels from distributing their noise throughout the compensated network. 

[0065] An other way to utilize this freedom is associated with the use of a reference subsystem of channels dedicated 
20 for compensation as in the prior art method 4. Although the method of this invention does not require such a subsystem 
it is still possible and in some situations useful to use one. In the present method these reference channels are treated 
in the same way as the actual measuring channels. In a case where the actual signal closely resembles an interference 
signal distribution in the measuring channel array (signal vector has a considerable projection on the interference signal 
space) it may be advantageous to use for the identification of the interference only the reference sensors which are 
25 immune to actual signals because they are typically located some distance away from the source of actual signals. On 
the other hand if there is a local disturbance seen by the reference channels only — noisy reference channel or vibration 
of the reference channel array, for example — the reference channels can be left out from the group of channels 
identifying the interference to prevent them from erroneously feeding the local disturbance into the measuring channels 
as a correction. 

30 [0066] The utilization of this freedom in the choice of cross coupling strengths is covered in claims 11 and 1 2. 

[0067] The actual cross-coupling matrix, which can be verified by experiment as described above, completely char- 
acterizes the state of the compensated multichannel device. It contains all the information about the system necessary 
for regeneration of the original uncompensated signals, if desired, from the measured data. Therefore, as regarding 
the data analysis, there are no complications even if the actual cross-couplings during the measurement would differ 

35 from the ones required for the ideal interference compensation. 

[0068] In summary, the approach of the invention enables to use sensitive magnetometer sensors for biomagnetic 
measurements instead of gradiometers, conventionally used to exclude environmental background interference. Al- 
though this particular embodiment of the invention has been described in detail for purposes of illustration, the same 
cross-coupling method can be used to make any type of A/-channel detector system immune to an external interference 

40 that is confined to a low dimensional subspace (n<N) of the signal space. 

Claims 

^5 1. A method for eliminating background interference from multichannel signal measurement, characterized in that 
outputs of at least two of the signal channels are cross coupled into a stage before the output of said signal channels 
in such a way that the background interference is minimized from the output of said channel. 

2. A method as defined in claim 1 , characterized in that 

50 

a compensating channel is arranged between the outputs of at least two of said signal channels and the stage 
before the output of said signal channel, 

the outputs of said signal channels are coupled to said compensating channel, and 

said compensating channel is coupled into the stage before the output of said signal channel. 

55 

3. A method as defined in claims 1 or 2, characterized in that the strength for coupling the output of said signal 
channel to said compensating channel is adjusted according to the calibration procedure. 
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4. A method as defined in claims 1-3, characterized in that the strength for coupling said compensating channel 
into the stage before the output of said signal channel is adjusted according to the calibration procedure. 

5. A method as defined in claims 1-4, characterized in that the calibration procedure comprises the steps of 

measuring the background interference by using the multichannel signal measurement, 
determining interference modes from said measurement, and 
- determining the strengths for said couplings on the basis of the interference. 

6. A method as defined in claims 1 - 5, characterized in that said strengths are determined by measuring the inter- 
ference without said cross coupling. 

7. A method as defined in claims 1-6, characterized in that said strengths are determined by measuring the inter- 
ference without test object. 

8. A method as defined in claims 1 - 7, characterized in that said strengths are determined by statistical analysis. 

9. A method as defined in claims 1-7, characterized in that said strengths are determined by principal component 
analysis. 

1 0. A method as defined in claims 1-7, characterized in that said strengths are determined by independent component 
analysis. 

11. A method as defined in claims 1-10, characterized in that said strengths are verified by exciting a single signal 
channel at a time with a spesific device. 

12. A method as defined in claims 1-11, characterized in that the number of said compensating channels equals the 
number of interference modes to be eliminated. 

13. A method as defined in claims 1-12, characterized in that said output channel is disconnected from said cross 
coupling if one of the predefined parameters of said channel has a value outside the predefined limits. 

14. A method as defined in claims 1 - 13, characterized in that the output of said signal channel is recorded as a 
function of time. 

15. A method as defined in claims 1 - 14, characterized in that the output of said compensating channel is recorded 
as a function of time 

16. A method as defined in claims 1-15, characterized in that the output of said signal channel is compared with the 
output of the compensating channel thus enhancing the output of said signal channel. 

17. A method as defined in claims 1-16, characterized in that the measured signal is a biomagnetic signal from the 
neuro- or cardiography source. 

18. A method as defined in claims 1-17, characterized in that said signal is measured by a SQUID-magnetometer. 

19. A multichannel signal measuring device, comprising a plurality of measuring channels provided with sensors, char- 
acterized in that the device comprises means for cross coupling outputs of at least two of the signal channels into 
a stage before the output of said signal channels in order to minimize background interference from the measure- 
ment. 

20. A measuring device as defined in claim 19, characterized in that said device comprises 

a compensating channel between the outputs of at least two of said channels and the stage before the output 
of the said channel 

a coupling element for coupling the output of said signal channel to said compensating channel, 
a coupling element for coupling said compensating channel to the stage before the output of said signal chan- 
nel. 
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21. A measuring device as defined in claims 19 or 20, characterized in that said coupling element cpmprises means 
for adjusting the strength for coupling the output of said signal channel into said compensating channel according 
to the calibration procedure. 

5 22. A measuring device as defined in claims 19-21, characterized in that said coupling element compnses means 
for adjusting the strength for coupling said compensating channel into the stage before the output of said signal 
channel according to the calibration procedure. 

23. A measuring device as defined in claims 19-22, characterized in that said device compnses means for 

10 

measuring the background interference by using the multichannel signal measurement, 

determining interference modes from said measurement, and 

determining the strengths for said couplings on the basis of the interference. 

15 24. A measuring device as defined in claims 1 9 - 23 characterized in that said device comprises means for determining 
said strengths by measuring the interference without said cross coupling. 

25. A measuring device as defined in claims 19-24, characterized in that said device comprises means for deter- 
mining said strengths by measuring the interference without test object. 

20 

26. A measuring device as defined in claims 19-25, characterized in that said device comprises means for deter- 
mining said strengths by statistical analysis. 

27. A measuring device as defined in claims 19-26, characterized in that said device comprises means for deter- 
25 mining said strengths by principal component analysis. 

28. A measuring device as defined in claims 19-27, characterized in that said device comprises means for deter- 
mining said strengths by independent component analysis. 

30 29. A measuring device as defined in claims 19-28, characterized in that said device comprises means for verifying 
said strengths by exciting a single channel at a time with a spesific device. 

30. A measuring device as defined in claims 19-29, characterized in that said device comprises equal number of 
said compensating channels and interference modes to be eliminated. 

35 

31. A measuring device as defined in claims 19-30, characterized in that said device comprises means for discon- 
necting said output signal channel from said cross coupling if one of the parameters of said signal channel has a 
value outside the predefined limits. 

40 32. A measuring device as defined in claims 19 - 31 , characterized in that said device comprises means for recording 
the output of said signal channel as a function of time. 

33. A measuring device as defined in claims 1 9 - 32, characterized in that said device comprises means for recording 
the output of said compensating channel as a function of time. 

45 

34. A measuring device as defined in claims 1 9 - 33, characterized in that said device comprises means for comparing 
the output of said signal channel with the output of the compensating channel thus enhancing the output of said 
signal channel. 

50 35. A measuring device as defined in claims 1 9 - 34, characterized in that the measured signal is a biomagnetic signal 
from the neuro- or cardiographic source. 

36. A measuring device as defined in claims 19-35, characterized in that said device comprises a SQUID-magne- 
tometer. 

55 

37. A multi-channel magnetometer for measuring biomagnetic signals provided with a set of separate channels for 
interference compensation, and applying the interference suppression method described in claims 1-18, char- 
acterized in that the elements providing the cross coupling between the sensor channels get their input as a 
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properly weighted linear combination from a selectable but arbitrary subset of the totality of all channels in order 
to compensate any dominant interference mode. 

38. A device described in claim 37, characterized in that said weights for each channel contributing to the input of 
5 the cross coupling elements are chosen in an optimal way based on a statistical analysis of the noise level of each 

channel. 

39. A method for collecting multichannel signals comprising of the signal of interest and of superposed background 
interference contributions which may be much larger than the signal of interest, characterized in that the mul- 

10 tichannel measuring device is provided with one or several elements that cross couple all or part of the signal 

channels before the readout of their output in such a way that the output of the channels in response to the back- 
ground interference can be made zero while the information concerning the background interference components 
is recorded by said elements. 

15 40. A method for cross coupling the signal channels of a multichannel device as defined in claim 39, characterized 
in that the coupling strengths between the channels and said elements are determined from a measurement of 
the background interference made with said multichannel device itself. 

41. A method as described in claims 39 and 40, characterized in that the cross coupling strengths between the chan- 
20 nels of the device are obtained by applying principal component analysis or independent component analysis, or 

any other method of statistical analysis on said measurement of the background interference. 

42. A multichannel device utilizing the cross coupling method described in claim 39, characterized in that the elements 
that record the background interference components and cross couple the signal channels before readout of the 

25 output are components in an analog electronic circuit, said registers obtaining input from and contributing to the 

output of the channels that they cross couple. 

43. A multichannel device utilizing the cross coupling method described in claim 39, characterized in that the elements 
that record the background interference components and cross couple the signal channels before readout of the 

30 output are registers in a digital signal processor, said registers obtaining input from and contributing to the output 

of the channels that they cross couple. 

44. A multichannel device utilizing the cross coupling method described in claims 39 and claims 42 or 43, character- 
ized in that the channels operate with the principle of negative feedback and said cross couplings between the 

35 channels are realized by exchange of negative feedback between the channels, said exchange of feedback being 

mediated by said cross coupling elements. 

45. A multichannel device according to claim 43, characterized in that the feedback providing the cross coupling is 
summed up with the ordinary negative feedback and delivered via the ordinary negative feedback path, or delivered 

40 separately via a separate feedback path. 

46. A multichannel device having mutually cross coupled signal channels according to claim 39, characterized in that 
the channels and cross coupling elements can be excited individually for the purpose of checking and re-adjusting 
said cross coupling strengths. 

45 

47. A multichannel device utilizing the cross coupling method described in claims 39 and 42, characterized in that 
the data on the state of the cross coupling elements as a function of time is collected simultaneously with the output 
data from the channels and stored on a memory device along with the output data. 

50 48. A device utilizing the method of claims 39 and 40, and the statistical analysis of the background measurement 
defined in claim 41 , characterized in that the number of elements activated to provide the cross couplings can be 
chosen on the basis of said statistical analysis. 

49. A device utilizing the method of claims 39 and 40, and the statistical analysis of the background measurement 
55 defined in claim 41 , characterized in that the elements activated to provide the cross couplings can be chosen to 

correspond to an arbitrary linear combination of the principal components of interference found in said statistical 
analysis. 
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50. A multi-channel magnetometer for measuring biomagnetic signals provided with a set of separate channels for 
interference compensation, and applying the interference suppression method described in claims 39 - 42, char- 
acterized in that the elements providing the cross coupling between the sensor channels get their input as a 
properly weighted linear combination from a selectable but arbitrary subset of the totality of all channels. 

51. A device described in claim 49, characterized in that said weights for each channel contributing to the input of 
the cross coupling elements are chosen in an optimal way based on a statistical analysis of the noise level of each 
channel. 

52. A method for measuring biomagnetic signals with a multi-channel magnetometer provided with a set of separate 
channels for interference compensation, and applying the interference suppression method described in claims 1 
- 18, characterized in that the elements providing the cross coupling between the sensor channels get their input 
as a properly weighted linear combination from a selectable but arbitrary subset of the totality of all channels in 
order to compensate any dominant interference mode. 

53. A method described in claim 52, characterized in that said weights for each channel contributing to the input of 
the cross coupling elements are chosen in an optimal way based on a statistical analysis of the noise level of each 
channel. 
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FIG. 2 



